Natural optical activity of metals 
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We derive the current response of noncentrosymmetric metal in normal and in superconducting 
states to the electromagnetic field with finite frequency and wave vector. The conductivity tensor is 
found contains the linear in wave vector off diagonal part generating natural optical activity. The 
Kerr rotation of polarization of light reflected from the metal surface is calculated. Its value is 
expressed through the fine-structure constant and the ratio of light frequency to the band splitting 
due to the spin-orbit interaction. 
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The natural optical activity or natural gyrotropy is 
well known phenomenon typical for the bodies having no 
centre of symmetry 1 . In this case, the tensor of dielectric 
permeability has linear terms in the expansion in powers 
of wave vector 



,(w,q) = eij(u,0) +ijijiqi, 



(1) 



where jn-l is an antisymmetric third rank tensor called 
the tensor of gyrotropy. The optical properties of a nat- 
urally active body resemble those of the magneto-active 
media having no time reversal symmetry. It exhibits dou- 
ble circular refraction, the Faraday and the Kerr effects. 

The description of the natural optical activity in terms 
of linear spacial dispersion of permeability 1 is appropriate 
for solid or liquid dielectric media. Whereas in the case 
of metals, it is more natural to formulate them in terms 
of spacial dispersion of conductivity tensor having the 
following form: 



;0,q) 



i(u,0) - iXijiqi- 



(2) 



The metals without inversion symmetry have recently 
become a subject of considerable interest arising mostly 
due to the discovery of superconductivity in CePtsSi. 2 
Now the list of noncentrosymmetric superconductors 
has grown to include UL£, CeRhSi 3 4 , CelrSisA Y 2 C 3 6 , 
Li2(Pdi_2;,Pt x )3l3i, KOs20e^, and other compounds. 
The theory has been mostly developed for description 
of superconducting properties of such a type of materials 
(for review see Ref. |9|). 

The spin-orbit coupling of electrons in noncentrosym- 
metric crystal lifts spin the degeneracy of the electron en- 
ergy band causing a noticeable band splitting. The Fermi 
surface splitting can be observed by the de Haas-van 
Alp hen effect discussed theoretically in the paper Ref. 
Hoi The band splitting reveals itself in the large residual 
value of the spin susceptibility of noncentrosymmetric su- 
perconductors at zero temperature^ 1 It also makes pos- 
sible the existence of nonuniform superconducting states 
those can be traced to the Lifshitz invariants in the free 
energy*^ 2 - Another significant manifestation of the band 
splitting is the natural optical activity. Here, we present 
the derivation of general expression for the current re- 
sponse to the electro-magnetic field with finite frequency 



and wave vector valid for the normal and the supercon- 
ducting state of noncentrosymmetric metals. We apply 
this result to calculation of Kerr rotation of polarization 
of light reflected from the surface of metal with cubic 
symmetry. In this case, the usual part of the conduc- 
tivity tensor is isotropic <7y (w, 0) = a(u>)Sij and the gy- 
rotropic conductivity tensor = Ae^/ is determined 
by the single complex coefficient A = A' + iX" such that 
a normal state density of current is 



o\E + A rotE. 



(3) 



We shall find here that the gyrotropy conductivity A is 
directly proportional to the ratio of the light frequency 
to the band splitting value. 

The current response to the electromagnetic field at 
finite q and u can be written following the textbook 
procedure 1 ^ generalized to two-band case ir>i^. In ap- 
plication to our situation, one has to remember that in 
a noncentrosymmetric crystal, all the values like single 
electron energy 



W(k) = £o(k)<5 a/ 3 +7(k)<T a/3 , 



(4) 



velocities v Q ^(k) = d£ Q/ 3 (k)/<9k, the inverse effec- 
tive mass (m~- 1 ) a p = <9 2 £ a/ 3(k) / 'dkidkj, the Green 

functions G q/3 (ti, k; r 2 , k') = -(T T ak a (n) a^fa)) and 
,F Q/3 (Ti,k;T2,k') = (T T a kQ (ri)a_ k /^(T2)) are matrices in 
the spin space (see for instance Ref. UA\ ). Taking this in 
mind, we obtain 



-e 2 7> 



IT,— — OO 



ii(w n ,q) 

d 3 k 



(5) 



{u i (k)G(°)( J Pr + )%(k)G(°)(JC_) 



(2tt) 3 



Aj(w n ,q). 



The transposed matrix of velocity is determined as 
v'(— k) = <9£*(— k)/<9k. The arguments of the 
zero field Green functions are denoted as K± = 
(f2 m ± u> n /2, k ± q/2). The Matsubara frequencies take 
the values Q m = ir(2m + 1 — n)T and to n — lixnT . We 
put here % = c = 1 and return back to the dimensional 
units in the final expressions. 
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Unlike to the paper |l6|, we do not write here the part 
of the current density arising due to the gauge shift in 
the argument of the order parameter A(k — eA(w n ,q)) 
describing magneto-optical phenomena in the supercon- 
ductors with broken time-reversal symmetry. We will 
return to this problem in a separate publication. 

Let us pass from the spin to the band representation, 
where the one-particle Hamiltonian 

H = ^^/3( k ) a ka a k/3 = ^( k ) C lcA C kA (6) 

k k,A=± 

has diagonal form. Here, the band energies are £ A (k) = 
£o(k) + A|"y(k)|. The diagonalization is made by the fol- 
lowing transformation 



(7) 



A=± 



with the coefficients 




V2I7KM+A7,)' 



forming a unitary matrix u(k). 

The zero field Green functions in the band representa- 
tion are diagonal and have the following form: 11 



where 



G^,(« n ,k)=J AA /G A (a; n) k), 
F^,(Lu n ,k)^S XX 'F x (uj n ,k), 

iuJ„ + £a 



(9) 



G A (w„,k) = - 



F\(u> n ,k) 



^ + ^ + |A A (k)p' 
i A (k)A A (k) 



and 



*A(k) 



;2+^ + |A A (k)P' 

7 a: (k) - ijyjk) 
7^(k)+7^(k)' 



(10) 



(11) 



The functions A A (k) are the gaps in the A-band quasipar- 
ticle spectrum in superconducting state. In the simplest 
model with BCS pairing interaction v g (k, k') = — V g , the 
gap functions are the same in both bands: A+(k) = 
A_(k) = A and we deal with pure singlet pairing 17 . 

Transforming the Green functions according the eqn. 
d?]) to the band representation, we obtain for the trace of 
matrices in eqn. ([5J): 

Tr{v t (k)G m {K+)Vj (k)G<°> (K-) 
+v i (k)FW{K + )v t j (-k)F+W(K_)} = 
v+ + Ak)G + v ++ j(k)G + +v ++li (k)F + v ++d (-k)Fl + 
w__,i(k)G_«__j(k)G_ + u__, i (k) J F_u__ J (-k)4 + 

w + _ )i (k)G_»_ +)i (k)G + + « + _ ii (k)f , _t;_ +)J -(-k)4 + 
V- +li (k)G + v + - d (k)G^+v- +ii (k)F + v + ^ j (-k)Fl.(l2) 



For the brevity, we have omit here the arguments of the 
Green functions. They are the same as in the upper two 
lines. The matrix velocity in the band representation is 



v AV (±k) = 



u 



a&(±k) 

9k 



, AQ (±k)v a/3 (±k) U/5A ,(±k) = 
y dk — Ti > AA '( k )' 



(13) 



where r(k) = u' (k)au(k) are hermitian matrices. We 
neglected the difference between u(k) and u(k ± q/2), 
which is of the order of (fcp<5) _1 , where 5 is the skin 
effect penetration depths 

For the calculation of gyrotropy of conductivity, only 
the last two lines in eqn. (fT2| consisting of interband 
terms are important. They are equal to 



dki dkj 



H+-T* +_[G_G+-F_M] 



r ; ,_ + r t ;,_ + [G + G_-F + Fl]}. 



(14) 



(8) Using the identity 



T l. + - T rn, + - = T l-+ T m-+ = <>lm ~ 77m + ieimnln , (15) 

where 7 = 7i/|"y|, one can rewrite the above expression 



dli dim 
dki dkj 

-F.Fh 4 



{(Sim ~ 77m)[G-G+ + G+G_ - (16) 



iei r 



il n [G-G+ - G+G_ - F-Fl + F + F f _}}. 



The linear part of current with respect to the wave 
vector q originates from the last term in the eqn. (|16[) . 

Of \ 9 9 ^ ^ f d k 

j,f(w„,q) = -te m e 7 T ^ J 7^3 x 



m— — 00 



*h[G-{K+)G-{KJ) - G+(K + )G-{K_) - 
F-(K+)Fl(K-) +F + (K + )Fl(K_)]A j (u> n ,q L ). (17) 

The momentum dependence of the spin-orbit coupling 
is determined by the crystal symmetry. For the cubic 
group G — O, which describes the point symmetry of 
Li2(Pdi_ x ,Pt a; )3B, the simplest form compatible with 
the symmetry requirements is 



7(k) = 7ok, 



(18) 



where 70 is a constant. For the tetragonal group G — 
C^v, which is relevant for CePtsSi, CeRhSi3 and CeIrSi3, 
the spin-orbit coupling is given by 

7(k) = l±(k v x - k x y) + "f\\k x k y k z (k 2 x - k^)z. (19) 

One can show that for the tetragonal group G — Gi„, 
the linear in the component of wave vector q part of con- 
ductivity is absent. So, we continue calculation for the 
metal with cubic symmetry where 7 = sigwjQ k. In the 
following, we put 7 = k taking 70 as a positive constant. 
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Let us find first the gyrotropy conductivity in the nor- 
mal state. Performing summation over the Matsubara 
frequencies, we obtain 



d 3 k 



kr x 



jf(w„,q) = -ieijie 7 j 

r /(e+(k_))-/(e-(k + )) 

l>„ + £+(k_)-£_(k+) 

/(£_(k_))-/(e f (k+)) 



e-(k_)-e+(k + ) 



Aj(o; n) q) 



(20) 



Here /(£±(k±)) is the Fermi distribution function and 
k± = k ± q/2. By changing the sign of momentum k — > 
k in the first term under integral and making use that 
£a (k) is even function of k, we come to 



j, 9 k,q) 



(2tt) 3 

/(e + (k + ))-/(e_(k_)) 

K)2 + (?+(k+) _ e _ (k _ )) : 



^^_(k_)-e+(k + )] 



r A 3 -(w„,q).(21) 



Analytical continuation of this expression from the dis- 
crete set of Matsubara frequencies into entire half-plane 
lj > is performed by the usual substitution iu) n — > LU+iS 
with <5 — ► 0. We see that the gyrotropic part of current 
is an odd function of the wave vector. At the same time, 
it is an even function of frequency. Hence, it presents a 
sort of displacement current originating of band splitting 
in noncentrosymmetric metal. 

We shall be interested in the low frequency and small q 
limit of this general formula. Expanding this expression 
up to the first order in q n u> 2 and performing integration 
over momentum space in the limit %u> <C 70 &f «Ef, we 
obtain 



J, ? (w,q) = -iei 



e 2 uj 2 



q n Aj(uj,q) 



(22) 



or, after substitution of the Fourier component of the 
vector potential by the Fourier component of an electric 
field A = E/iu, 



3i(u, q) = - e ijn 10 2 , q n Ej(u, q). (23) 



This corresponds to 



A 



127t 2 7oA:f 



(24) 



To calculate the Kerr rotation, we consider linearly 
polarized light normally incident from vacuum to the 
boundary of a medium with complex index of refraction 



N = n + in 



(25) 



expressed through the diagonal part of complex conduc- 
tivity a = a' + ia" by means of the usual relations 



1 



47TO-" 



2nn 



4-Ka' 



to 



(26) 



The light is reflected as elliptically polarized with the 
major axis rotated relative to the incident polarization 
by an amount^ 3 - 



(1 - n 2 + k 2 )Ak + 2nnAn 
= (1 - n 2 + k 2 ) 2 + (2nn) 2 ' 

where for the current given by eqn. ([3]). 

An = n + — ri— = , 



An 



4ttA' 



(27) 



(28) 



(29) 



are the differences in the real and imaginary parts of the 
refraction indices of circularly polarized lights with the 
opposite polarization. 

After substitution the eqn. (J24J) in eqns. (|29|). (|28| we 
find that Ak = and An expresses through ratio of the 
light frequency to the band splitting 2^okp as 



a 



Tilo 



3tt 7ofc.F 



(30) 



Here, a — e 2 /Tic is the fine structure constant. 

It is worth to be noted that the coefficient A consists of 
two parts: contribution of the quasiparticlcs in between 
the Fermi surfaces of two splitting bands and in the vicin- 
ity of the two Fermi surfaces (compare with Ref. [TTh . In 
the superconducting state at fiui ^> A, A can be given as 



A 



2 



87r 2 7ofci? 



1 



1 



Y(T) 



and 



HT) = /i 



2T cosh 2 ^ 2 + A 2 /2T) 



(31) 



d£ (32) 



is the Yosida function. At T — 0, the second term in eqn. 
(|32|) . contribution of the quasiparticles in the vicinity of 
the Fermi surfaces, disappears and only the first term 
remains. As a result, we obtain in superconducting state 
at T = 



An 



2ir 7ofci? 



(33) 



In conclusion, we presented here the derivation of the 
current response to the electromagnetic field with finite 
frequency and wave vector in noncentrosymmetric metal. 
The conductivity tensor contains a gyrotropic part re- 
sponsible for the natural optical activity. As an example 
the Kerr rotation for the polarized light reflected from 
the surface of noncentrosymmetric metal with cubic sym- 
metry is calculated. The found value of the Kerr angle 
is expressed through the fine structure constant and the 
ratio of the light frequency to the spin-orbit band split- 
ting. The result can be used for the direct experimental 
determination of the latter value. 
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